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Although crucial for resolving the issue of 
charge symmetry in the nuclear force, 
direct measurement of nn-scattering by 
coUiding free neutrons has never been per- 
formed. At present the Russian pulsed 
reactor YAGUAR is the best neutron 
source for performing sucb a measure- 
ment. It has a through channel where the 
neutron moderator is installed. The neu- 
trons are counted by a neutron detector 
located 12 m from the reactor In prelimi- 
nary experiments an instantaneous value 
of 1.1 X 10^^/cm^s was obtained for the 
thermal neutron flux density. The experi- 
ment will be performed by the DIANNA 
Collaboration as International Science & 
Technology Center (ISTC) project No. 
2286. 
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1. Introduction 

After many years of conflicting reports for the neu- 
tron-neutron ^Sq scattering length a^, by the late 1990s 
the results appeared to have converged to consistent 
numbers: a^ = -18.55 ± 0.05 (stat) ± 0.3 (syst) fin [1] 



from the ^H(7i;", Yii)n reaction, and ann = -18.7 =t 0.3 
(stat) ± 0.6 (syst) fm [2] from the ^H(n, nn)p reaction. 
Comparison with the well established value of the pro- 
ton-proton scattering length a^ = -17.3 ± 0.005 (stat) ± 
0.4 (syst) fm indicates charge- symmetry breaking 
(CSB). 
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However, a recent investigation of the neutron-neu- 
tron and neutron-proton final-state interaction in the 
n-d breakup reaction revived this issue. The a^n-value 
of -16.27 ±0.40 fin was obtained in Ref [3] for the 
^H(n, np)n reaction, where the total uncertainty consists 
of about equal parts of statistical and systematic errors. 
This latest result differs fiom that of Ref. [2] by almost 
4 (J, and renews the long-standing controversy regard- 
ing the conclusion on charge-symmetry breaking. The 
possible solution is a direct measurement of the nn- 
scattering length. Such a measurement has been never 
performed. Here we describe specifics of the direct nn- 
scattering measurement which is under preparation at 
the unique aperiodic pulsed reactor YAGUAR. 



2. General Scheme of the Experiment 

In effective-range theory, the singlet scattering cross 
section a^ is defined by the scattering length a^ as 



a = 47ra^ at ^ ^ 0, 



(1) 



where k is the neutron wave number. For thermal neu- 
trons the zero effective-range approximation works 
well. The effective thermal neutron cross section a 
measured with unpolarized neutrons is a statistical sum 
of the singlet o; and triplet a^ cross sections 



13 1 2 

a = —G^ +—G, =—G^ =7ia„n 

4^4*4^ nn 



(2) 



Since the Pauli exclusion principle for identical parti- 
cles forbids the interaction of two neutrons in the triplet 
state, the cross section a^ is expected to be zero, and the 
measured cross section a is equal to na^,. 

Proposals for a direct measurement of the nn-scatter- 
ing length have a long history; detailed references are 
provided in Ref [4]. None of these proposals has been 
implemented. Except for the proposed use of under- 
ground nuclear explosions, most of the proposals are 
similar. The neutrons collide in a neutron interaction 
chamber, which we denote as the nn-cavity, and the 
scattered neutrons are detected externally In such 
arrangements the "target" and the "beam" are neutrons 
produced by the same neutron source. Therefore, the nn 
scattering intensity is proportional to the square of the 
neutron flux density, while the background intensity 
depends linearly on the flux density. 

The situation is most favorable at pulsed neutron 
sources, where the thermal neutron flux density can be 
very high. It appears that the pulsed aperiodic reactor 
YAGUAR [5] at Snezhinsk, Russia, fulfills require- 



ments for the nn cross section measurement. This nn 
project was initially proposed at the Dubna ISINN-VIII 
meeting [6]. During the fast neutron pulse of 0.68 ms 
(FWHM) duration, about 10^^ fast neutrons are generat- 
ed in the active core. The 40 L volume of the YAGUAR 
active core contains the water soluted salt UO2SO4. The 
solution contains 465 g/L of uranium which is enriched 
to 90 % ^"^^U. The effective diameter of the cylindrical 
active core is 40 cm, with an effective inside diameter 
of 17 cm. The critical height is about 39 cm, depending 
on the position of the startup rods. The body of the 
reactor has a central channel of 1 5 cm diameter, which 
in the standard operating mode has startup rods that 
leave a cylindrical space of 12 cm diameter. 

A direct measurement of the nn-scattering length is 
to be performed by the time-of-flight (TOP) method by 
counting the scattered neutrons arriving after the reac- 
tor burst at a detector placed behind special collimators 
far away fiom the reactor. If only the colliding neutrons 
contribute to the detector counting rate, and if the 
parameters of the neutron field are determined and the 
geometry is known, then the detector counts (integrat- 
ed over the thermal neutron part of the TOP spectrum), 
measure the nn-scattering cross section. The proposed 
arrangement for the nn-scattering experiment at 
YAGUAR is shown schematically in Pig. 1 . The active 
core 3 is placed on three supports at 2 m above the floor 
level. The neutron CH2 moderator 4 is inserted inside 
the channel. An evacuated tube contains a collimation 
system 5, and the neutron detector 6 is placed at a flight 
path of about 12 m. The collimation system is designed 
to screen the source and to minimize background due to 
neutrons scattered fiom the walls. The "^He absorber 1 
reduces neutron scattering fiom the back wall. 



3. Neutron Flux 

The reactor produces a hard neutron spectrum in the 
through-channel (average neutron energy 0.91 MeV) 
that is not suitable for the nn-experiment. We per- 
formed calculations for various moderators in the 
YAGUAR geometry and also performed measurements 
with various moderators at the reactor YAGUAR. The 
Monte Carlo modeling [7] was performed using the 
code MCNP-4 [8]. The cylinder height of 40 cm was 
fixed and the thickness of the polyethylene cylindrical 
shell varied. 

Neutron activation studies were performed to meas- 
ure the thermal neutron densities inside cylindrical con- 
verters inserted into the YAGUAR reactor channel. The 
polyethylene converters were fabricated as hollow 
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Fig. 1. Setup for the nn-scattering experiment. The horizontal scale 
is ten times greater than the vertical scale. See text for explanation. 



cylinders of 39.7 cm length and 12.0 cm outer diame- 
ter, with different inner diameters. The activation detec- 
tors for the absolute flux measurements were gold and 
copper foils placed in the central plane. One detector 
set was covered with a 1 mm thick Cd screen. Details 
of the measurements are given in [4]. The relative dis- 
tribution of the neutron flux along the axis of the 
through-channel was also obtained. These measure- 
ments, combined with known neutron activation cross 
sections, determined the fluency, defined as the number 
of neutrons per cm^ per YAGUAR pulse. 

The measured neutron fluencies for the four polyeth- 
ylene converters vary linearly with thickness. This 



dependence agrees very well with Monte Carlo simula- 
tions [7]. With the experimental results for the fluence 
and using the value A^ = 0.7 ms for the thermal neutron 
pulse, we conclude that the 3 cm thick cylindrical mod- 
erator provides, for the pulse power of 30 MJ, an 
instantaneous value of 1.1 x lO^Vcm^s for the thermal 
neutron flux cP^entrai in the central region of the channel. 
The moderated neutron flux spectrum is expected to 
be predominantly Maxwellian and to have a 1/^" 
epithermal tail with the value of n depending on the 
moderator thickness. The gold and copper activation 
data were analyzed to obtain experimental estimates of 
such parameters as the effective Maxwellian tempera- 
ture, the relative size of the epithermal tail, and the 
slope parameter n. The neutron flux spectrum generat- 
ed with these parameters is shown in Fig. 2 for the 2.2 
cm thick moderator. The MCNP modeling of the neu- 
tron fluxes for different moderator thicknesses agrees 
with the experimental activation results. Based on these 
results we chose the moderator thickness value of 2.5 
cm for future nn measurements. 
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Fig. 2. Neutron flux F{E) versus neutron energy E inside the cylin- 
drical cavity with the 2.2 cm thick moderator. The spectrum is nor- 
malized to unity at the maximum. 



4. Simulations 

4.1 Counting Rates 

After nn-scattering in the evacuated moderator chan- 
nel, some of the scattered neutrons will reach a detec- 
tor placed at about 1 2 m from the channel central plane, 
and will be measured by the TOF method with an effec- 
tive solid angle 13^^. For ideal collimation, direct neu- 
tron paths from the moderator and the walls (besides 
the back wall of the nn cavity) to the detector are 
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excluded. Then the detector counts A^^n P^r pulse (inte- 
grated over the thermal part of the TOP spectrum) are 
related to the nn-scattering cross section a^, the aver- 
age neutron flux density 0^^, the effective pulse dura- 
tion At, the effective nn-cavity volume V, and the most 
probable velocity Vq by 

^nn = 2c,, ^a^AtVq^ (counts/pulse), (3) 

where the constant c^v is determined by Monte Carlo 
and analytical calculations for a given moderator geom- 
etry and neutron flux distribution. For a pure thermal 
neutron spectrum and including the CM/LAB transfor- 
mation of the scattered flux in the direction of the 
detector, the value of c^^ = 0.83 ± 0.01 [9]. Background 
considerations and the effects of fast-thermal and fast- 
fast collisions are discussed in subsections that follow. 
With a realistic value for the solid angle (Q^^ - 5 x 
10"^), then from Eq. (3) we expect N^ ~ 180 nn-counts 
in the neutron detector per burst. 

4.2 Backgrounds 

Of paramount importance for the direct nn-scattering 
measurement is the detector background. Ideally with 
extensive shielding no slow or fast source neutrons can 
reach the detector directly. However, in practice one 
can expect several background components. 

The thermal neutron background component from 
the thermal and epithermal neutrons multiply scattered 
on collimators will depend on the collimation system. 
Several collimator options, as presented in [4] and in 
references therein, were considered, and an estimate of 
this component was found to be less than 20 % even for 
the simplest two-diaphragm collimation. With the use 
of Boron- 10 coated collimators one can reduce this 
component considerably. The inevitable scattering 
from the back wall will be reduced by choosing an opti- 
mal geometry and the use of the TOP technique. 

Fast neutrons that originate during the reactor pulse 
are also separated by the TOP measurement. The time 
interval allocated for detecting thermal neutrons is 
about 5 ms and is delayed by ^2 ms relative to the max- 
imum of the reactor pulse. At this time the subcritical 
reactor continues to emit fission neutrons. During this 5 
ms interval about A'^- 10^^ neutrons from the bottom 
area of the reactor vessel are directed into the lower 
hemisphere. These neutrons penetrate partially through 
shields and can present a background problem. This 
issue was addressed in Ref [10]. The Monte Carlo code 
MCNP4-C with the weight window variance reduction 
technique was chosen to calculate the neutron transport 



for a 6 m model of the shield consisting of water, con- 
crete, and earth, with a 6 cm diameter channel lined 
with boron carbide. Scaled to the real geometry the 
results obtained correspond to a fast neutron back- 
ground level of ~5 %. 

The background-to-nn signal ratio from neutron 
scattering on the residual gas depends on the gas pres- 
sure and does not depend on the collimation system. At 
a gas pressure of <10"^ torr, this background component 
is estimated to be <1 %. 

Therefore, preliminary calculations on the back- 
ground issues are promising, and more extensive mod- 
eling by transport code experts is underway at 
Snezhinsk. Direct experimental separation of the nn- 
signal from the background is also possible, because 
the nn-scattering intensity is proportional to the square 
of the neutron flux density, while the background inten- 
sity depends linearly on the flux density. 

4.3 Time-of-Flight Spectra 

The nn scattering measurements at the pulsed reactor 
YAGUAR will be performed by the TOP method. The 
transformation from the energy spectrum of the inci- 
dent flux F(E)dE (as shown in Pig. 2) to the TOP spec- 
trum A^(^) is accomplished by the use of equations 
F(E)dE = N(t)dt and dE = -2(dt/t)E. Por a 12 m flight 
path and 2.5 cm thick cylindrical moderator the simu- 
lated "incident" TOP spectrum is shown in Pig. 3 by the 
full curve. Such a spectrum would be measured by the 
neutron detector if it were placed as shown in Pig. 1 and 
the nn cavity filled with a heavy gas such as argon, for 
which the transformation of the neutron spectrum after 
scattering can be neglected. The spectrum essentially 
consists of a Maxwellian part above ^ = 2 ms and a l/f' 
epithermal part. There is also an intermediate distribu- 
tion for t between 1.6 ms and 2.0 ms. Por convenience 
we refer to the distributions above and below ^* = 2 ms 
as the "slow" and "fasf components. This value of ^* 
is expected to be a cut off value for the experimental 
TOP spectrum below which the background (not shown 
in Pig. 3) will increase sharply. 

However, for the collision of light particles with 
equal masses, energy and momentum conservation 
leads to the well known "step" distribution of final 
energies extending down to zero and to an anisotropy 
of scattering in the laboratory (LAB) coordinate sys- 
tem. Questions arise about the shape of the scattered 
neutron spectrum in the YAGUAR experiment and con- 
cerning the size g of a possible contribution to the ther- 
mal energy region from "fast-slow", g|, and "fast-fasf , 
g^, collisions producing "slow" neutrons. These ques- 
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Fig. 3. Calculated time-of-flight spectra N{t) (in relative units) for the 12 m flight path. Full line: the "incident" 
(before scattering) neutron flux from the 2.5 cm thick cylindrical moderator; broken line: the nn-scattered flux. 
The spectra are normalized to equal areas. 



tions can be answered in two ways: analytically by 
invoking a complicated double differential scattering 
cross section in the LAB system, or numerically by 
calculating the collision rates for initial neutron trajec- 
tories (for neutrons emitted from the internal cylindri- 
cal surface of the moderator), simulating the isotropic 
angular distribution of scattering in the CM coordinate 
system, and finally by performing analytical transfor- 
mation of the scattered velocity vectors to the LAB sys- 
tem in order to follow neutron trajectories to the detec- 
tor. 

For the latter modeling, a new code was written as 
described in [11]. The TOF spectrum of the scattered 
neutrons calculated with this code is shown in Fig. 3 by 
a broken curve while the initial non-scattered spectrum 
is drawn by the full curve. With these spectra normal- 
ized to equal areas the difference in shapes is quite evi- 
dent. However, the thermal parts of the scattered spec- 
tra for cases of the realistic and the Maxwellian inputs 
remain nearly the same. The code performs a proper 
analysis of the "fast-fast", "fast-slow" and "slow-slow" 
contributions to the thermal part of the realistic spec- 
trum. According to this analysis, the integral of the real- 
istic scattered spectrum above ^* changes only slightly 
as compared to the pure Maxwellian case. The prelim- 



inary numbers are: gfs~6%, ^<0A%; graphical 
illustrations are given by Crawford [11]. 

These issues are relevant to the c^^ value in Eq. (3) 
which is now expected to change only slightly between 
the cases of the pure Maxwellian and realistic spectra. 
Calculations of c^ for the YAGUAR neutron flux spec- 
trum with an epithermal tail are currently in progress. 
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